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1 Introduction 

Event-related desynchronization (ERD) is 
attenuation in the amplitude of oscillations in 
specific frequency that occurs in relation to a 
specific neural activity [1]. The opposite 
phenomenon, event-related synchronization (ERS), 
is an increase in the amplitude [2]. A previous MEG 
study demonstrated that ERD occurs in alpha and 
beta bands following electrical stimulation of the 
median nerve or self-paced finger movements in 
human [3]. However, the changes in brain rhythm in 
higher frequency bands following somatosensory 
stimulation have been reported only in the animal 
brain [4]. Furthermore, the volumetric distribution 
of those changes has not been demonstrated 
previously. 

In the present study, we have developed a novel 
analysis method, synthetic aperture magnetometry 
(SAM) [5], to reveal the frequency-dependent 
volumetric distribution of somatosensory evoked 
magnetic fields (SEFs) on the individual MR images. 
SAM is one of the statistical spatial filtering 
methods using an adaptive beamformer, analogous 
to the technique which obtains the high spatial 
selectivity from radio antenna arrays [5,6,7]. 

2 Materials and Methods 

2.1 Data acquisition 

Seven healthy right-handed subjects (six males and 
one female, age 27-35 years) participated in this 
study. Informed consent for the study was obtained 
from all subjects. During the experiment, the 
subjects sat in a comfortable chair in a magnetically 
shielded room with eyes open. The right median 
nerve was electrically stimulated with a constant 
voltage, square-wave pulse of 0.2 ms duration 
delivered transcutaneously at the wrist. Stimulus 
intensity was adjusted to the lowest level to produce 
a twitch of the thumb. The interstimulus intervals 
varied randomly from 2 s to 4 s. SEFs were recorded 
with a 64-channel whole-head MEG system with 


third-order SQUID gradiometers (Whole-Head 
Squid System Model 100, CTF Systems Inc., Port 
Coquitlam, B.C., Canada) [8]. MEG signals were 
digitized at a sampling rate of 1250 Hz and filtered 
with a 400 Hz on-line low pass filter. Notch filters 
were used at 60 Hz and 120 Hz to eliminate AC line 
noise. The data of 1200 ms duration with a 600 ms 
pre-stimuli interval were collected for each of 100 
trials. The head position was registered with the 
localization coils at the nasion and bilateral 
preauricular points at the beginning and end of each 
session. For each subject, magnetic resonance 
imaging (MRI) was obtained with an 1.0 T MRI 
system (MAGNETON impact, SIEMENS Inc., 
Japan) in a T1-weighted sequence of 130 sagittal 
slices (1.5 mm thickness). The nasion and bilateral 
preauricular points were identified on MR images 
with the aid of fiducial markers, allowing the 
overlay of MEG data on the individual MR images. 

2.2 SAM analysis 

The MEG data were filtered with each of the 
following passbands: alpha (8-13 Hz), beta (13-25 
Hz), low gamma (25-50 Hz) and high gamma (50- 
200 Hz). The region of interest (ROI) was set to 
include the whole brain with a 5 mm voxel 
resolution. A spherical area of 1.5 cm diameter 
located at the center of head frame was excluded 
from the ROI because the signal noise increases at 
the area remote from the SQUID sensor. Then the 
signal power of each voxel was estimated by SAM 
[5]. Student’s t values were calculated for the 
changes in the signal power in each voxel between 
the active states after the stimulus and the control 
states before the stimulus. Time windows were 
selected as follows; 1) 10 to 100 ms vs. -100 to -10 
ms, 2) 10 to 600 ms vs. -600 to -10 ms, and 3) 100 
to 600 ms vs. -600 to -100 ms. Only 25-50 Hz and 
50-200 Hz bands were analyzed for the latency 10- 
100 ms taking account of the Nyqwuist’s sampling 
theorem. The distribution of Student's t values was 
displayed on the individual MR images. 
Localization of the spectral changes in signal power 



before and after stimulus was accessed in the voxels 
where Student's t value indicates the maximum 
value. 

The power spectra of signal sources in the voxels 
corresponding to the hand area of the bilateral 
primary somatosensory cortices (S-I) were 
compared between before (-600 to -10 ms) and after 
(10 to 600 ms) the stimulus. 

2.3 Dipole modeling 

Data were also analyzed with conventional single 
equivalent current dipole (ECD) modeling in order 
to assess the consistency between SAM analysis and 
ECD analysis. After eliminating the artifact- 
contaminated trials, the data were passed through a 
150 Hz off-line low pass filter and averaged. The 
pretrigger 100 ms were used for DC offset. The 
location of dipoles was calculated according to a 
least-square fitting algorithm for 600 ms starting at 
the onset of the trigger [8]. The dipoles with the 
least fitting errors were chosen, excluding those with 
a fitting error higher than 30%. 


3 Results 

3.1 SAM analysis 

In all subjects, prominent signal power increase, or 
ERS, was demonstrated in the high gamma band 
(50-200 Hz), which was localized strictly around the 
contralateral central sulcus (FIG 1A). The maximum 
t value in the time window of 10-600 ms was 
10.69±2.73. In addition, ERD was also observed in 
the alpha (8-13 Hz), beta (13-25 Hz), and low 
gamma bands (25-50 Hz) bilaterally around the 
central sulcus in 4 subjects, (FIG 1C, D). In one 
subject, ERD was demonstrated only in the alpha 
band (8-13 Hz) around the contralateral central 
sulcus. In the remaining 2 subjects, ERD was 
demonstrated only in the low gamma band (25-50 
Hz) around the ipsilateral central sulcus. 

The ERS in 50-200 Hz was observed mainly within 
100 ms after the stimulus (FIG 1A). No major ERS 
in 50-200 Hz was demonstrated later than 100 ms 
after the stimulus (FIG IB). However, ERD in 8-13 
Hz, 13-25 Hz, and 25-50 Hz was sustained later than 
100 ms after the stimulus. 



FIG 1: SAM statistical images of a representative subject showing frequency-dependent volumetric 
distribution of Student's t value between the signal powers before and after the transcutaneous electrical 
stimulation of the right median nerve. The event-related synchronization (ERS) in high gamma band (50-200 
Hz) was localized strictly around the contralateral primary sensorimotor area in 10-100 ms (A). No major 
ERS in 50-200 Hz was demonstrated in 100-600 ms (B). The event-related desynchronization (ERD) in the 
alpha (8-13 Hz) (C), beta (13-25 Hz) (D) and low gamma bands (25-50 Hz) were distributed bilaterally 
around the primary sensorimotor area whereas the latter also involved the posterior parietal area. The 
equivalent current dipole (ECD) corresponding to N20m located at the contralateral primary somatosensory 
area is indicated in red (A). 
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FIG 2: The power spectrum of signal source power 
in the hand area of the right (A) and left (B) primary 
somatosensory cortex before (-600 to -10 ms) and 
after (10 to 600 ms) the stimuli. Before the stimuli, 
the power spectrum revealed biphasic peaks in the 
alpha and beta bands (black line). After the stimulus, 
both peaks were attenuated demonstrating the ERD 
in these bands (red line). 

ERD in 8-50 Hz in the posterior parietal area was 
also demonstrated in 4 subjects (FIG ID), and ERS 
in 8-13 Hz in the left premotor area in 2 subjects 
(FIG IB). In 4 subjects, either ERD or ERS was also 
localized around the upper bank of the Sylvian 
fissure. 

The power spectra of signal sources in the voxels 
corresponding to the bilateral SI had biphasic peaks 
in the alpha and the beta bands before the stimulus. 
Both those peaks were attenuated after the stimulus, 
corresponding with the bilateral ERD on the 
individual MR image (FIG 2). The changes in the 
power spectrum in high gamma band were not 
apparent between the control and active states. 
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FIG 3: The waveforms of the somatosensory evoked 
magnetic fields. Within 80 ms after the stimulus, 
typical early components, N20m, P27m, P38m, and 
P51m were consistently identified (arrowheads). 

3.2 Dipole modeling 

Within 80 ms poststimulus, typical early 
components (N20m, P27m, P38m, and P51m) were 
identified consistently in the filtered and averaged 
MEG waveforms (FIG 3). The ECD for each of 
those components, N20m, P27m, P38m, and P51m, 
was localized in the posterior wall of the 


contralateral central sulcus, corresponding to the 
hand area of the contralateral S-I area (FIG 1A for 
N20m). The average fitting error for ECDs of N20m 
was 13.10±9.98 %. ECDs were inconsistently 
localized around the upper bank of the Sylvian 
fissure later than 80 ms poststimulus. 

3.3 Consistency between the location of 

ECD and ERS 

The deviation of the location of the ECD N20m 
from that of ERS in high gamma band (50-200 Hz) 
in the time window of 10-600 ms was 4.6-15.7mm 
(7.4 ±3.8 mm). 

4 Discussion 

The present study demonstrated the significant 
signal power attenuation (ERD) in the alpha and 
beta bands around the bilateral central sulcus 
following the stimulation of the right median nerve. 
The characteristics of this ERD agree well with the 
previous results obtained with EEG and MEG 
studies [3,9]. They are localized bilaterally in SM-I 
with predominance to the contralateral hemisphere, 
with main frequency in the alpha and beta bands, 
occurred immediately after stimulation and 
sustained for more than several hundreds of 
milliseconds [3]. 

The SAM statistical imaging also demonstrated the 
prominent ERS in the high gamma band (50-200 
Hz), strictly localized around the contralateral 
central sulcus mainly within 100 ms after the 
stimulus. However, this change was not clear with 
the power spectrum analysis. The failure with the 
power spectrum analysis might be related to the high 
background noise in this frequency band. Several 
studies have also demonstrated the high frequency 
responses in the human sensorimotor system. The 
EEG study by Pfurtscheller et al demonstrated 40 
Hz ERS over the contralateral central region during 
self-paced finger movement [10]. They interpreted it 
as related to the process of motor planning. Curio et 
al. and Hashimoto et al. revealed the extremely high 
frequency activities (600 Hz) in the contralateral S-I 
after the electrical stimulation of the median nerve 
[11,12]. They supposed that those activities were 
evoked by the GABAergic inhibitory intemeurons in 
area 3b [12]. Many animal studies reported the 
gamma band oscillation around 40 Hz in the 
somatosensory system [4,13] as well as in the other 
sensory systems such as olfactory and visual 
systems [14,15]. The present study revealed the ERS 
in the human primary somatosensory system in the 
high gamma band (50-200 Hz) previously not 





























reported. 

The waveforms of the MEG data contained four 
main components within the first 80 ms, N20m, 
P27m, P38m, and P51m, whose ECDs were 
localized in the posterior wall of the contralateral 
central sulcus. A number of investigators suppose 
that N20m is representative of the activation of 
pyramidal cells in area 3b, whereas the sources of 
the later components are still controversial 
[16,17,18]. The S-I [16] primary motor area [17], 
and other sensory or motor areas [18] are the 
possible generator of these components. Since these 
components form the frequency higher than 50 Hz, 
it is possible that they are involved in the ERS in 50- 
200 Hz demonstrated in this study. It is also possible 
that a part of this high frequency ERS can be 
attributable to the activation of pyramidal cells in 
the primary sensorimotor area. 

Both ERS and ERD are closely related to the brain 
signal processing on specified brain activation. In 
the present study, frequency-dependent volumetric 
analysis demonstrated the significant difference in 
the spatiotemporal distribution of ERS and ERD in 
different frequency bands: ERDs in 8-13 Hz, 13-25 
Hz, 25-50 Hz were distributed bilaterally around the 
SM-I, also tended to involve the posterior parietal 
area, and sustained for more than 100 ms. In contrast, 
ERS in 50-200 Hz was localized strictly in the 
contralateral SM-I and attenuated within 100 ms. 
ERS in 8-13 Hz and 13-25 Hz was observed in the 
premotor area on the contralateral side. This 
significant difference in the spatiotemporal 
distributions of ERS and ERD may support the 
hypothesis that the neural mechanism for ERS and 
ERD are different from each other, reflecting the 
responses of different cell assemblies related to the 
specific neural functions rather than the frequency 
shift of the same cell assembly. 

Conclusion 

A novel MEG analysis (synthetic aperture magneto- 
metry, SAM) revealed a frequency-dependent 
volumetric distribution of SEF after electrical 
stimulation of the right median nerve. The 
prominent feature was a short-lasting high- 
frequency ERS in 50-200 Hz strictly localized on 
the contralateral primary sensorimotor cortex 
whereas ERD in lower frequency bands (8-13 Hz, 
13-25 Hz, and 25-50 Hz) was observed bilaterally 
and sustained later than 100 ms. This is the first 
report that demonstrated the difference in the 
spatiotemporal distribution of ERS and ERD as well 
as the difference in the frequency bands of them in 
the human primary sensorimotor system following 


peripheral nerve stimulation. 
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